Radiative atom-atom interactions in optically dense media: Quantum corrections to 

the Lorentz-Lorenz formula 

Michael Fleischhauer and Susanne F. Yelin 
Sektion Physik, Ludwig-Maximilians Universitdt Miinchen, D-80333 Miinchen, Germany 
and: Dept. of Physics, Texas A & M University, College Station TX 77843-4242, U.S.A. 

(February 1, 2008) 



Generalized single-atom Maxwell-Bloch equations for op- 
tically dense media are derived taking into account non- 
cooperative radiative atom-atom interactions. Applying a 
Gaussian approximation and formally eliminating the degrees 
of freedom of the quantized radiation field and of all but a 
probe atom leads to an effective time-evolution operator for 
the probe atom. The mean coherent amplitude of the local 
field seen by the atom is shown to be given by the classical 
Lorentz-Lorenz relation. The second-order correlations of the 
field lead to terms that describe relaxation or pump processes 
and level shifts due to multiple scattering or reabsorption of 
spontaneously emitted photons. In the Markov limit a non- 
linear and nonlocal single-atom density matrix equation is 
derived. To illustrate the effects of the quantum corrections 
we discuss amplified spontaneous emission and radiation trap- 
ping in a dense ensemble of initially inverted two-level atoms 
and the effects of radiative interactions on intrinsic optical 
bistability in coherently driven systems. 



I. INTRODUCTION 

The interaction of the radiation field with a dilute en- 
semble of atoms is usually described in the semiclassical 
and dipole approximation by the well-known Maxwell- 
Bloch equations. This description fails to be accurate, 
however, when a dense medium is considered. 

Since the early work of H. A. Lorentz and L. Lorenz |l]] 
it is known that the classical local field, that couples to 
an atom in a dense medium, differs from the macroscopic 
(Maxwell) field by a term proportional to the medium 
polarization ^|J|]. The most prominent effects of the 
Lorentz-Lorenz (LL)-correction in dense media are the 
change of the linear index of refraction according to the 
Clausius-Mossotti relation , the enhancement of non- 
linear susceptibilities$s M , shifts and deformation of res- 
onance lines dJ(|, intrinsic optical bistability fflM, and 
piezo-photonic switching |J. 

On the other hand the quantum nature of the ra- 
diative atom-atom interaction can drastically influence 
the behavior of the ensemble. In the extreme case of 
anisotropic, high-density samples, excited atoms can co- 
operatively emit spontaneous photons, a phenomenon 
known as superradiance |pT5|— p~2| . But even if the system 
does not fulfill the conditions for cooperative evolution, 
the presence of spontaneous photons and the associated 



effects like amplified spontaneous emission (or superlu- 
minescence) and radiation trapping Jl3[ can not be ne- 
glected. Imprisonment of incoherent photons especially 
affects otherwise long-lived ground-state coherences. We 
therefore expect radiative atom-atom interactions to be 
important in areas such as resonant linear and nonlinear 
optics based on atomic phase coherence |lj,|l^], cooling 
of atoms and Bose-Einstein condensation via velocity- 
selective coherent population trapping |fl6|| and optical 
computing. 

Another important effect of large atomic densities is 
the increase of atomic collisions. Here we will not con- 
sider these effects, however, and focus our attention en- 
tirely on radiative interactions. 

In the present paper we study the atomic evolution in a 
dense medium irradiated by external coherent light fields. 
The macroscopic classical radiation field in the medium 
obeys Maxwell's equations with the mean atomic polar- 
ization as source term. To derive equations of motion 
for the many-atom system, we start from a nonrelativis- 
tic quantized interaction Hamiltonian. Thus interactions 
between the atoms mediated by the quantized radiation 
field such as reabsorption and scattering of spontaneous 
photons are taken into account. 

Our aim is to derive an effective single-atom density- 
matrix equation. For this we introduce an interaction 
picture with the radiation field coupling to all other 
atoms. Assuming a Gaussian (and therefore classical) 
statistics of the interacting field, we can formally elim- 
inate its degrees of freedom from the probe-atom time 
evolution. In the Markov limit of short-lived field cor- 
relations this yields a density-matrix equation for the 
probe atom. We will show that the mean coherent am- 
plitude seen by the probe atom differs from the macro- 
scopic Maxwell field by a term proportional to the mean 
polarization of the medium in agreement with the clas- 
sical Lorentz-Lorenz relation [0. In addition, the den- 
sity matrix equation contains relaxation and level-shift 
terms, which describe reabsorbing and multiple scatter- 
ing of spontaneously emitted photons. The correspond- 
ing relaxation rates and frequency shifts are proportional 
to the spectrum of the incoherent part of the radiation 
inside the medium. This spectrum is also the Fourier- 
transform of a certain 2-time Greensfunction, for which 
we derive a Dyson equation. A formal solution of the 
Dyson equation allows to express the incoherent spec- 
trum in terms of atomic variables. Thus we eventually 
obtain a closed, nonlinear and spatially nonlocal density 
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matrix equation of Lindblad-type. 

Our paper is organized as follows. In Sec. II we de- 
rive the effective single-particle time-evolution operator 
by formally eliminating the degrees of freedom of the 
quantized radiation field interacting with the background 
atoms. In Gaussian approximation this operator contains 
first and second-order field cumulants. In the Markov 
limit of spectrally broad field correlations, a density ma- 
trix equation is obtained. In Sec. Ill we show that the 
first-order term leads to the Lorentz-Lorenz relation be- 
tween the coherent amplitude of the local field, the mean 
field amplitude in the medium (Maxwell field), and the 
mean polarization. In Sec. IV we derive a Dyson equation 
for the second-order field cumulants or 2-point Greens- 
functions and formally solve them in terms of single-atom 
density matrix elements. The resulting nonlinear density 
matrix equation is discussed in Sec.V for the examples of 
amplified spontaneous emission and radiation trapping 
in an inhomogeneously broadened system of initially ex- 
cited two-level atoms and intrinsic optical bistability in 
a strongly driven dense ensemble of two-level atoms. 



II. EFFECTIVE TIME-EVOLUTION OF ATOMS 

A. Formal elimination of the quantized radiation 
field 

We here consider an ensemble of atoms interact- 
ing with the quantized radiation field under conditions 
which justify the dipole and rotating-wave approxima- 
tion (RWA). Since we are interested in the dynamics of 
a single atom, we distinguish a probe atom at position 
r*o with a dipole operator p and environment atoms at 
positions Tj whose dipole operators are denoted by d? . 
The Hamiltonian of the system is given by 



H = J2 H o+ Afield - P ■ \E(r ) + £{f ) 



Y^d 1 -[E{r 3 )+£(r 

j¥0 



(1) 



where Hq and -fffioid are the free Hamiltonians of the jth 
atom and the quantized radiation field respectively, and 
we have split the field in an operator component E and an 
external classical driving field £ . We use an interaction 
picture where the time evolution is described by 

S = Tev{-tfjrV 9 {T)} = 

Texp | ~J dr p(r) [E(f , r) + £ (f , r)]| , (2) 

were T denotes time-ordering and the field operator E 
still contains the coupling to all other atoms. For nota- 
tional simplicity we have suppressed vector indices of the 
dipole moment and electric field. With the help of (0) any 



(time-ordered) correlation function of probe-atom oper- 
ators Ah and Bh in the Heisenberg-picture (subscript 
" H" ) can be related to interaction picture operators via 



{T- l [A H {t x )A H {t 2 )]T[B H {H)B H {U)\) = 
(T- 1 [5- 1 A(t 1 )^(t 2 )]T[5 J B(t 3 )S(t 4 )]), 



(3) 



where (• • •) stands for Tr{po ■ ■ •} with po = p{— oo) being 
the initial density operator at t — —oo. 

A very helpful formal simplification of Eq.(|J) can be 
achieved by introducing the so-called Schwinger-Keldysh 
time contour C JTtJ shown in Fig. 1 which starts at 
t = — oo, goes to t — +oo and back t — —oo. Each 
physical time correspond two times on the contour. A 
time ordering operator Tq is introduced, which is iden- 
tical to T on the upper branch (+) and to T^ 1 on the 
lower branch (— ) of the contour and orders all operators 
with time arguments on (— ) to the left of those with time 
arguments on (+). 



+ 
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FIG. 1. Schwinger-Keldysh time contour 

With these definitions we can write Eq. (||) with a single 
exponential time-evolution operator. This will consider- 
ably simplify the following elimination procedure. 

{[T^AH^AHitiWBnit^BHiU)]) = (4) 
(Tc[S c A(ti)Afe)B(tt)B(ti)}), 

where the superscripts ± specify the branch of the con- 
tour, and 



Sc 



T, 



c exp 



df V p (f) 



(5) 



with f denoting a time on C . 

We now formally eliminate the degrees of freedom of 
the quantized radiation field and the environment atoms 
by tracing over the corresponding states. In order to 
express the expectation value of an exponential operator 
again as an exponential operator, i.e. as a new - effective 
- time-evolution operator, we use a generalization of the 
cumulant generating function for a classical stochastic 
variable X filpl: 



exp{ s A}) =exp £ ^((A" 1 )) 



(G) 



where the ((X m )) are the cumulants, which have the fol- 
lowing explicit form 



((XY)) = (XY) - (X)(Y), etc. 



(7) 
(8) 
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As can bee seen from (g) , the elimination procedure leads 
in general to an infinite number of terms in the effective 
action. To make the problem tractable, we will however 
assume that the radiation field is Gaussian, i.e. that all 
cumulants ((E m )} with m > 2 vanish identically. This is 
a consistent and for our purposes well justified approxi- 
mation. With this we find 



S C = ( s c) fie i d 



T c expj - jdrp{f) £(f ,f) + (E(r ,f)) 



j df i J df 2 p(f 1 )T>(f ,fi;r ,f 2 )p(f 2 ) 



(9) 



where 



2V(1,2) = Ut c £^(ri,n)K(r 2 ,f 2 ; 



(10) 



is a (tensorial) Greensfunction (GF) of the interact- 
ing electric field, and we have used the abbreviations 
f = fx,Ti and 2 = f 2 ,r 2 . Note, that we used a 
short notation, and p(l)T>(i,2)p(2) in Eq.(|) should 

read Y^,,p=i P<*(t) VaffO-, 2) pp(i, 2). We now apply the 
rotating-wave approximation. For this we introduce 
slowly varying positive and negative frequency compo- 
nents, 

p(f) = p +(f) +p-(f) = p + {f)e- WT + p-(f)e luJT , (11) 
E(f) = E+(f) + E-(f) = E+{f)e- lUT + E-(f)e luJT , (12) 



B. Markov approximation and single-atom density 
matrix equation 

The effective single-atom time-evolution operator dl3| ) 
leads in general to integro-differential equations of mo- 
tion. We therefore restrict the discussion to situations 
that justify a Markov approximation, i.e. we assume 
that the characteristic decay time of field cumulants is 
short compared to the characteristic time of the atomic 
dynamics. This is the case, for example, if the medium 
is inhomogeneously broadened or under quasi-stationary 
conditions. We note that the Markov approximation used 
in the present paper rules out superradiance, since the su- 
perradiance time is shorter than the decay time of field 
correlations [ [il"||l2j ]. In order to describe fast cooperative 
processes, pair-correlations between different atoms need 
to be taken into account as a dynamical variable. This 
will be discussed in a future publication, where we derive 
an effective density matrix equation for atom pairs [ pcfl . 

The Markov approximation amounts to assuming a 8- 
correlation of D^ v and C^ t „ in physical times. 



D^(t,t')=D^(t)5(t-t') 



(17) 
(18) 



with A, B e {+, — } explicitly denoting the contour 
branches. It is convenient to introduce dimensionlcss 
dipole operators c, , such that p+ (t) = p p (t) and 
p~(t) ~ p^cr^i) (and corresponding relations for the 
slowly- varying variables) . With this we eventually arrive 
at 



with lo being the transition frequency of the considered 
probe atom, and neglect combinations of the type p + E + 
and p~ E~ . Thus we have 



oett _ 



T c exp< - Jdf p + {f)£ L (r ,f) +p (f)££(r ,f) 



I dfi I df; 



„.a / ^'l / "-'2 

2h 2 Jc Jc 



P + (n) D(r Q ,f 1 ;r ,f 2 )p (f 2 ) 



+p (f 1 )C(r ,f 1 ;ro,f 2 )p + (f 2 )]|, (13) 

where 

£ L ^t) = £^(r,t) + (E^t)) (14) 
is the local field seen by the probe atom, and 

((TcE-ifo^E+ifo^e-^-^, (15) 
C^(r , fi;r ,f 2 ) = 

((T c E+(r ,n)E-(^f 2 ))) e+Mn-r,). (lg) 



S c c s T c exp< 



dr <J ti {T + )£ Lll {rQ,T) - a f ,{T^)£ Lfl {f ,T) 

Xj 

+al(T + )£+ fi (f Q ,T) - ai(T-)£^(r Q ,T) 



0>( r -) £ 4( T "-) - 2cr Al (r_)crt( T+ ) 



(19) 



00 / r^(cj,r) J^(u,t) 
-oo V 2 2 
<4(t_)o- m (t_) - 2o-t( T _) - At ( T+ ) 



4( T +)°>( T +) + 



+ hj dTH t^( UJ ,T) |ct m (r+ ) al (r+ ) - 

0" i uO~-)<4(' r -) - 4(t+)ct a1 (t+) + at (r_)o- M (r_) 



+ - / dTh^(uj,T) ct+(t + )(J m (t + ) - o-t(T_)o>(Y_) 



The lower indices ± at the time argument denote the 
branch on the Schwinger-Keldysh contour which is rele- 
vant for operator ordering under the action of Tq- The 
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first term in (|19|) describes the interaction of the probe 
atom with the local field in RWA. 



([E;(f ,t),E+(f ,t-T)})( 
([E-(r ,t),E+(r ,t + T)}}e 1 ' 



r 

^ / dr((E-(f ,t)E+(f ,t + T))), 



(25) 



(20) 



is a positive hermitian matrix, whose eigenvalues describe 
decay and pump rates induced by the incoherent pho- 
tons inside the medium. Eq. (|20|) has a simple physi- 
cal interpretation. The incoherent radiation inside the 
medium causes stimulated transitions from excited to 
ground states and vice versa. The corresponding rate 
is proportional to the spectral density of the radiation 
taken at the atomic transition frequency. Apart from 
some dimensional constants D~J~ is precisely the spec- 
tral energy density of the incoherent field at the position 
fo and at the transition frequency u> of the probe atom. 

lfll/ (u, t) = ^ j^dr([E+(r , t + r), E~ (f 0) *)]) - + ^ | 



is the corresponding spontaneous contribution. Within 
the approximations made, reflects the Lamb-shift of 
excited states altered by the presence of the medium. 
Here we are not interested in the Lamb shift and therefore 
consider it included in the free Hamiltonian H . 

The effective time-evolution operator ([l9|) directly 
leads to the following master equation for the single-atom 
density operator: 



^l^ + ^l^P 



(26) 



~Y + ~Y ) \ vIv/mP + P<A,°~ii - 2o- M pcrJ 



(21) 

is the spontaneous contribution to the "down rate" in 
the atomic medium. (Note that the commutator con- 
tains the field operators interacting with the environment 
atoms.) Since we are not interested here in the effects of 
the medium to the spontaneous decay, we replace 7^„ 
by the free-space value 7°^. We will show in Appendix 
A, that Eq.(|2l]) indeed leads to the well-known Wigner- 
Weisskopf result for radiative decay in free space, if we 
replace E by the free field. Light-shifts induced by the 
incoherent component of the radiation field inside the 
medium are described by the hermitian matrix 



This is the first main result of the present paper. We note 
that this equation is nonlinear and nonlocal, since the 
light-shift and decay matrices depend via the field cor- 
relations on the surrounding atoms. The equation does 
however have the Lindblad form pl| and thus preserves 
positivity and the total probability. In order to obtain 
a closed set of equations, we calculate in the following 
sections the yet unknown quantities £l , r M „ , and in 
terms of single- atom density matrix elements. 



H^ u (uj,t) 
h 2 



dr 



((E-(f ,t)E+(f ,t-r)))e- 



III. THE AVERAGE LOCAL FIELD AND THE 
LORENTZ-LORENZ RELATION 

We recognize from Eq.(p6|) that the probe atom is cou- 
pled to a classical (c-number) field of amplitude 



((E-(f Qi t)E+(f ,t + T)))e +iur 



Eq.(p2|) can also be expressed in terms of D + : 



(22) 



£ L (f,t)=£(f,t) + (E(f,t)). 



(27) 



2nh 



, D-+(f ,uj';t) 



2tt 



du' 



UJ — U) 

T ILU {u',t) 



id 



(23) 



where P denotes the principle part of the integral. In 
systems with inhomogeneous broadening the collective 
light-shifts are often negligible as they are usually small 
compared to the inhomogeneous width. 



The first term is the external coherent field ( = field in the 
absence of the medium), and the second term is the mean 
coherent amplitude of the field scattered by all other 
atoms. Note, that the contribution of the probe atom it- 
self is not included. On the other hand, the macroscopic 
field £m, which enters Maxwell's equations is the total 
field inside the medium (averaged over a spatial region 
large compared to the characteristic atomic distance, but 
smaller than A 3 ). Thus the local field, given in ( p7| ) dif- 
fers from the macroscopic Maxwell field essentially by 
the scattering contribution of the probe atom itself. In a 
continuum approximation we find 



h 2 



dr 



(24) 



£La(r,t) = £ Ma (r,t)~ 

% -q [ d 3 f' f°°dt' Dl%(f,t;r',t'){j> Hff (t')), 

n JK C J-00 



(28) 
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where Dq ' is the free-field retarded propagator and (jph) 
is the expectation value of the probe-dipole operator (in 
the Heisenberg picture), g is the atomic density and K t 
denotes integration over a small sphere of radius e. The 
retarded propagator of the electric field is given by p2fl 



(29) 



if i 



9(r) 



dx^dx^ 



S(r — ct) 



Iri — r*2| and O is the Heaviside 



where r = ti 
step function. 

When substituting _Dg et f rom Eq.(p9|) into Eq.@) we 
note that in the limit e — > only the term which results 
from the second spatial derivative of 1/r survives. Using 



d 2 



dx$dx$ •' 



~ = -y ^ (3) ( 



ft - fa) <W 



(30) 



we find 



£ L {r, t) = £ M (f, i) + g^fi(fe (*) 
= 5m(F,<) + — P 



(31) 



which is identical to the classical Lorentz-Lorenz rela- 
tion E10 when we identify V = q(ph)- It should be 
mentioned that the Lorentz-Lorenz relation holds for the 
mean amplitude of the field and not for the field opera- 
tors itself as claimed in pjj . 

Making use of (|3l]) we can define an effective semiclas- 
sical interaction operator 

VL = -Y,i%(t)i L(t (f j ,t) . (32) 



IV. QUANTUM CORRECTIONS 

We now discuss the light-shift and decay matrices in 
the generalized density-matrix equation ( |26| ) in more de- 
tail. Both depend on the field cumulants or Greensfunc- 
tions 



D-+(r,t;r,t') = ({E-(r,t)E+(r,t'))). 



(33) 



(Note that the superscript " — h" indicates that the first 
time argument is on the lower and the second time ar- 
gument on the upper branch of the Keldysh contour and 
has nothing to do with the frequency components of the 
field.) The aim of the present section is to calculate D h 
in terms of atomic variables. For this we apply non- 
equilibrium Greensfunction techniques [p4[. 



1. Dyson-equation for D(l, 2) 

We define the exact and the (known) free Greensfunc- 
tions (GF) on the Keldysh contour as 

ZV(1,2) = «T c £;(f 1 ,t 1 )S+(f 2 ,t 2 ))), (34) 
D ^(l,2) = ((TcE^ftME+fai*))), (35) 

where E$ denotes the free field, i.e. without coupling 
to the medium. The contour-Greensfunction D(l,2) 
contains four real-time GFs: D++(l,2), D~ + (l,2), 
D + -(l,2), and D — (1,2), where the superscripts "±" 
specify contour branches. The first and the last are the 
time- and anti-time ordered propagators and the retarded 
and advanced propagators are given by the combinations 



£> rot (l,2) =U ++ (1,2)-D + -(1,2) 
= D-+(l,2)-D—(l,2), 

D adv (l, 2) = D++{1, 2) - D-+{1, 2) 
= £>+-(!, 2) -£>— (1,2). 



(36) 
(37) 



Within the KWA and in the absence of thermal photons, 
we have 



D++ (1, 2) « 0^(1,2), 
0^(1.2) «0, 

D+- p (l, 2) « D^(l, 2) - Dltp(l, 2), 
^(1,2)^-^(1,2). 



(38) 
(39) 
(40) 
(41) 



A formal solution to the atom-field interaction can be 
given in terms of a Dyson- integral equation p4[ , by in- 
troducing a formal polarization function Tl a p(l, 2) 

ZV(1,2) = D ^(l,2)- 

JJjtV d2'D^ a {\, V) n^(l', 2') Df3 V (2', 2). (42) 

Here f c dl denotes integration over the Schwinger- 
Keldysh contour as well as spatial integration over the 
medium. The Dyson-equation ( |42"| ) represents nothing 
else than a formal summation of the perturbation se- 
ries where the polarization function is determined by the 
medium response. We now have to find a good approxi- 
mation for 11(1, 2). 

2. self-consistent Hartree approximation 

One easily verifies that in lowest order in the atom-field 
coupling, the polarization function is given by a correla- 
tion function of dipole operators of non- interacting atoms 



(43) 



I0U) = 

£ f?E<( T ^(*0^(*2)» £ree 5(f 1 -f^Cfa -P S ). 
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This corresponds to a Hartree approximation in many- 
body theory. This approximation is physically justified, 
when the nonlinear light-shift and decay terms do not af- 
fect the atomic dynamics, that is if the probability that a 
specific atom reabsorbs or scatters a spontaneous photon 
is small. Such a situation is realized, for example, in the 
classical case of radiation trapping where a small number 
of photons (much smaller than necessary to saturate the 
medium) is trapped in a dense absorbing medium [ p^[ . 
We are here also interested, however, in situations, where 
incoherent photons significantly alter the atomic dynam- 
ics. A consistent approximation, which accounts also for 
these cases is the self-consistent Hartree approximation, 
where the cumulants of free dipole operators in (43) are 
replaced by cumulants of interacting dipole operators. 

n Q/3 (i,2) = 

^ E<< T c°ia(*ito/J&)>) rj) S(f 2 - r 3 ). (44) 

j 

As shown in Appendix B the Dyson-equation for the 
contour GF can be rewritten in the RWA in terms of the 
real-time GFs as follows: 



d3d4D^(l,3)n^(3,4)D$?(4,2), 

(45) 



where D™£(1, 2) (= Df^(2, 1)) obeys the Dyson equa- 
tion 

D%(1,2) = IC(1 S 2)- 



d3 dA DI%(1, 3) IC(3, 4) D$(4, 2). (46) 



Here the time integration goes from -oo to oo and we 
have introduced the atomic source correlation 



KAri,h;r 2 ,t 2 ) 



£«4.(*i)M*2)» S(n - r,) 6(r 2 - f 3 .) (47) 



as well as the atomic response function 



'-e(h-t 2 ) 



n 



The names reflect the physical meaning of the terms. 
The Fourier-transform of n s is proportional to the spon- 
taneous emission spectrum of the atoms and that of IF et 
gives the susceptibility of the medium. 



& 



D" + (l,2) 
D^d,2) 
D ret (l,2) 

n s (l,2) 



n ret (i,2) 



/ d r T 

FIG. 2. Feynman diagrams, definitions 




FIG. 3. Graphical representation of Eq.(^). The incoher- 
ent intensity at the position of the probe atom is the sum of all 
spontaneous contribution propagated through the medium. 



FIG. 4. Graphical representation of Dyson equation (^) 
for retarded GF inside the medium. Iteration generates 
all-order scattering contributions. 



Eqs.([45[) and (46) can be given an instructive graphical 
representation shown in Figs. 2-4. Eq.(|45j) (illustrated in 
Fig. 3) says that the incoherent radiation intensity is ob- 
tained by summing the spontaneous-emission contribu- 
tions from all atoms propagated through the medium. 
The iteration of the Dyson equation ( |4^ ) (shown in 
Fig. 4) describes multiple scattering of spontaneous pho- 
tons by atoms during the propagation from a source atom 
to the probe atom. 



3. Explicit expressions for the collective decay rate and 
light- shift 

We now approximately solve the Dyson-equation ( |4"6| ) 
for the retarded propagator in the medium. We first 
introduce a continuum approximation. 



d A r P™(r, h,t 2 ) S(n - r) 8{r 2 - r), 



(49) 



(50) 
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where N is the atom density and the overline denotes av- 
eraging over some inhomogeneous distribution. Similarly 



nL(fi,ti;f 2 ,t 2 ) 



d 3 r P*(r, h,t 2 ) S(fi - f) S(f 2 - r), 



(51) 



P^MM) = ((al( tl )a jv (t 2 ))) (52) 

Thus Eq.(||) reads 

DTp{?iM;r2M) = Dl e *p(n,ti;?2,t 2 )- (53) 

dt[ / dt' 2 / d 3 f i Dg*, (rl , ; f i , tj ) 

OO J — OO J V 

To solve this integral equation, we now make the fol- 
lowing approximations. We first extend the spatial in- 
tegration to infinity, which basically means that we are 
solving for the retarded propagator in an infinitely ex- 
tended medium. Secondly, we replace r' x in the atomic 
response function by r 2 , i.e. we evaluate the response 
at the position of the source. We furthermore consider 
quasi-stationary conditions, i.e. assume that P ret (ti,i 2 ) 
depends only on the time difference r = t[ — t' 2 - We 
only keep an overall slow (parametric) time dependence. 
This means that we consider propagation times short 
compared to the characteristic time of the atomic evo- 
lution, which is consistent with the earlier Markov ap- 
proximation. With these simplifications we can turn the 
integral equation ( |53| ) into an algebraic one by Fourier- 
transformation with respect to x = r*i — f 2 and r = t\—t 2 . 
At this point a word of caution is needed: As will be dis- 
cussed in Appendix C, the retarded GF in an amplifying 
medium is not Fourier-transformable, since it grows ex- 
ponentially with r — \x\. We therefore should view the 
transformations as finite-time and finite-space Fourier- 
transforms, and hence the algebraic equation as an ap- 
proximation. 

Using the definition 



D^ix^-t) 



ihu e q o r 
67reoc 2 r 



-iq a r 



(56) 



where A is the wavelength of the transition under consid- 
eration in the rest frame, and r = \x\ = Iri — f 2 \. 



Qo = q' Q {r, w, t) + iq'^f, w, t) = - 



l + ^P ret (r,u;;*) 
3e 



(57) 



q'n is the inverse absorption/amplification length in the 
medium and q' characterizes the corresponding phase 
shift. We here have assumed that |lm[P ret ]| 7i/3eo < 1. 

With Eq.(p6|) we can now express D h (r ,cj;<) in 
terms of atomic variables 



D-+(?o,«,t) 



fc2 4 



(6^c 4 J v 



d 3 f 



s 2q' ( l(r,u;t)r 



■P s (r,Lu;t). 
(58) 



Here r = \ f— fo\ is the distance between source and probe 
atom. With Eq.(|58|) we finally find for the collective de- 
cay rate and light-shift 



r(w,t) 

H(u,t) 



P 2 ^ 



(6^) 2 ^ 4 Jv 



e 2q'd (r,u;t)r „ 

d 3 f - P s (r>;i), (59) 



(6^ c 4 J y 



d 3 i 



dw' eVo'C^'itJr ps(jf )W;t ) 



2tt 



(60) 



Eqs.(p9|) and ( |60[ ) are the second major result of the 
present paper. In applying these results to a specific 
problem, we still have to calculate the source-correlation 
P s in terms of density matrix elements. This then yields 
a closed nonlinear and nonlocal density matrix equation. 
We will illustrate this for some examples in the following 
section. 



F(&«)= / d 6 x I dT¥{x,T)e~ lUT e lq ' x , (54) 

J Voo J —OO 

the solution of Eq. (f53|) reads 



D (q,tu;t) 



1+D (g»-P re >;i) 



D (q,u). 
(55) 



Here D ret and II rot denote 3x3 matrices in coordinate 
space and 1 is the unity matrix. 

For simplicity we now disregard polarization, i.e. we 
replace the 3x3 matrices by simple functions. We note 
however, that a generalization is straight forward. As 
shown in detail in Appendix C, we eventually arrive at 



V. EXAMPLES 

A. Inhomogeneously broadened two-level system 

We here consider an inhomogeneously broadened dense 
ensemble of randomly polarized two-level atoms in a 
cylindrical geometry as shown in Fig. 5. For this system 
the time-evolution of the dipole operator a = \b)(a\ is 
determined by the simple Heisenberg-Langevin equation 



&j = -(ico 3 ab + T ab )aj + noise, 



(61) 



where the noise term denotes a white noise source, which 
is however of no interest here. Lo° ab — ui® b + Aj is the 
atomic transition frequency in the laboratory frame. We 
here take into account Doppler-broadening which leads 
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to a shift Aj of the lab-frame transition frequency from 
the rest-frame frequency uP ab . The collective light-shift is 
small compared to the average Doppler-shift and there- 
fore neglected. The coherence decay rate r a f, consists of 
two contributions, one resulting from free-space sponta- 
neous decay 7 and the other from the collective decay T, 
Tab = T + 7/2. Eq.(|6l|) can easily be solved by Laplace- 
transformation (x(s,t) :— j dr e~ ST x(t + r)), which 
yields 



«at( S ;i)^-(i)» 



PlJt) 



S ~ iuJ ab + T ab ' 
3 



s + iui 3 ab + r ab 



From this we immediately obtain 



h T ab + i(oj -w£ 6 ) 



5s/-. n 2 P 2 », pia(t)T ab 



h' (r afc )2 + (c-a;i b )2 



(62) 
(63) 

(64) 
(65) 



where p is the dipole moment of the transition and the 
overbar denotes averaging over the velocity distribution 
of the atoms, which is given by the Gaussian distribution 



W(A 3 ) 




(66) 



FIG. 5. Dense sample of inhomogeneously broadened 
two-level atoms in cylindrical geometry. 

Since the lab-frame atomic transition frequency de- 
pends on the velocity, the collective decay rate, which 
is proportional to the incoherent radiation spectrum at 
this frequency, will be velocity dependent as well. Thus 
we have in general a set of nonlinear coupled equations 
corresponding to different velocity classes. If there are 
fast velocity-changing collision, the population dynamics 
of all velocity classes will however be approximately the 
same. In this case we may set 



not hold if the inhomogeneous broadening mechanism is 
not due to Doppler as for example in solids. In that 
case one has to consider the full set of equations. Using 
Eq. ( |67|) , we find in the limit of large Doppler- broadening 
Ad > r ab 



P ret (f,u,t) 



where A 



P 2n pr Paa(r,t) - pbb(r,t) 



'2tt 



-A 2 /2Af 



l-UI- 

7T 



A D 
2 £VA. 



dye 



-v 2 /2 



, (68) 



o — Lo® b is the detuning from the atomic 
resonance at rest. For the collective decay rate only the 
real part of P lct is important which enters the absorption 
coefficient according to Eq.([57|) 



p 2 N uj 



(69) 



3Ueo c 



— Paa{r,t) - Pbb(r,t) _A 2 /2A 2 n 



Similarly we have for the source term in the strong 
Doppler-limit 



P B (r,u,t) 



2p 2 N 



27T P*^,t) e _A 2 /2A? 



(70) 



Combining (j39|) and QTy) and applying the relation be- 
tween the free-space radiative decay rate 7 and the dipole 



moment p: 



37rhe c s ^f/uj A [E5f (cf. also Appendix A) 



yields the collective decay rate fl59| ) for a probe atom with 
(lab-frame) transition frequency u> at position r*o 



r(w,t) = 

7 / d 3 f2q%(r,u>,t) 



(71) 



e2q'o(r,u>,t)r 



Paa(r,t) 



47TT- 2 Paa(r,t) - p bb (r,t)' 

with r = \f — r*o|. To obtain the effective decay/pump 
rate we have to averaged over the velocity distribution 



r(i) 



duj ■ 

—00 



1 

/2vrAn 



(72) 



We now discuss two limiting cases. In the first case 
we assume a small excitation in the medium. This cor- 
responds to the classical situation of radiation trapping 
in an inhomogeneously broadened two-level medium. We 
will show that in this case Eq.(f72|) leads to the integral 
equation of Holstein |l3| . In the second case we will disre- 
gard the spatial dependence but keep the nonlinearities, 
and consider the temporal evolution from an initially ex- 
cited ensemble. 



(67) 



where fj, € {a, b} and r denotes the position of the atoms 
considered. Note, however, that this approximation does 



1. linear limit and Holstein equations of radiation trapping 

For small excitation, the retarded light propagation 
can be regarded as propagation in a medium with all 
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population in the lower state, i.e. p bb = 1 and p aa = 0. 
Thus 



n ,z ,s in ,n P^NuV^tt . 

g o(r ,U;,t)=g ( W) i) = - — e 



-A 2 /2A 2 D 



(73) 



and we can approximate the denominator in ( |71| ) by —1. 
This results in 



r(t)w 7 d 3 r G(r ,r)p aa (r,t) 
Jv 



where 
G(f ,f) 



(74) 



(75) 



2<Z V, t) 



2q"(u,t)r 



Aitr 2 



1 



4 7rr 2 



1 



47rr 2 / <9r 



(9 



— exp 



-Koe 



Here ifo = NX 2 g and 5 = j/V^nAu characterizes the 
ratio of the homogeneous to the inhomogeneous width. 

The dynamical evolution of the ensemble is described 
by the Bloch equation 



Paa{r Q ,t)=r(t)- 7 + 2r(i) p aa {r ,t) 



(76) 



In the small-excitation limit, the term Tp aa is of second 
order and can be neglected. We thus arrive at the linear 
integral equation for the atomic excitation 



paa(r ,t) 

-7Paa(r6,t)+7 I d 6 r G(r a ,r) p aa (r,t). 

IV 



(77) 



(78) 



Eq.(f77|) is the integro-differential equation for radiation 
trapping derived by Holstein in jl3| for the special case of 
Doppler-broadened two-level atoms. Thus in the linear 
limit we have rederived the theory of radiation trapping 
of @. 



2. dynamics of initially inverted two-level system in 
small-sample approximation 

Let us now discuss a nonlinear problem, but in a small 
volume, such that the space dependence can be disre- 
garded. In this case we can carry out the volume integral 
placing the probe atom on the axis of the long cylindrical 
sample (see Fig. 5). We find for the decay rate for a probe 
atom with transition frequency lj 



r(w,t) 



Paa(t) 



Pbb(t) ~ Paa(t) 



1 — exp 



(-K(t)c 



-A 2 /2A 2 D 



(79) 



where A = u> — u>® b , and 

K(t) = K d\p bb (t) - p aa (t)]. 



(80) 



Averaging over the inhomogeneous velocity distribu- 
tion of the atoms yields 
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Paa{t) 



Pbb(t) 



-j= I dye v ' 



Paa{t) 

1 — exp 



-K(t)e 



(81) 



Note that T(t)/j remains finite at p aa = 1/2, since the 
diverging denominator is multiplied by a vanishing inte- 
gral expression. 

The time evolution of the excited-state population 
from an initially completely inverted system is shown in 
Fig. 6 for different values of the density rj = NX 2 d =100 
(solid line) and r\ = 500 (dashed line) and g = 0.01. 
The dotted line corresponds to the free-space decay. One 
recognizes a non-exponential behavior, with an acceler- 
ated decay in the initial phase corresponding to amplified 
spontaneous emission and a substantial slow-down of de- 
cay in the long-time limit. 



R„ 



Tl= 

100 
500 



Y t 

FIG. 6. Time dependence of excitation in spatially homo- 
geneous, dense two-level medium. Time is in units of the 
inverse free-space decay rate, n = N\ 2 d = (dotted), 100 
(line), 500 (dashed), g = j/V^ttAd = 0.01. 

The effective rate of decay of the excitation r e g = 
—paa/ paa is shown in Fig. 7. One can see that for rj = 500 
the initial decay rate is already of the order of the in- 
homogeneous Doppler- width (logA^/7 ps 1.6) and the 
Markov-approximation of slow atomic evolution becomes 
invalid. For higher atomic densities the system would 
show superradiant decay in the initial phase which can- 
not be described by the single-atom density matrix equa- 
tion. As noted before, modeling of the cooperative decay 
requires a two-atom density matrix description, which 
will be discussed elsewhere [E0l. 
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yt 

FIG. 7. Effective rate of decay r e fi = —paa/Paa for ex- 
amples of Fig. 6. Dotted line corresponds to free-space de- 
cay. Amplified spontaneous emission in the initial phase and 
slow-down of decay in final phase are apparent. 

One also verifies from Fig. 7 for the case rj = 100, that 
the decay becomes exponential again in the long-time 
limit. The asymptotic escape rate is given by 



7esc — 



7 



K (n In Ko)V 2 ' 



(82) 



which can be orders of magnitude smaller than 7. This 
result agrees with Eq.(l.l) of Jl^]b up to a numerical fac- 
tor of the order of unity, which is due to the fact that we 
here have disregarded a possible spatial inhomogeneity. 

It is also instructive to consider the time-dependent 
spectrum of incoherent radiation or equivalently T(uj,t). 
This is done in Fig. 8 for 77 = 500. Shown is the spectral 
distribution at different times normalized to the averaged 
rate T(t). 



rtco.t) 




co-a)° b /A D 

FIG. 8. Spectral distribution of incoherent decay as func- 
tion of time. r/ = 500, g = 0.01. Dotted line shows 
Doppler-distribution of atomic transition frequencies (not 
normalized). Time is in units of 7 -1 . 

The dotted line shows the (not normalized) inhomo- 
geneous distribution of atomic frequencies according to 
Eq. (|66|) . One recognizes that the incoherent spectrum 
broadens with the decay of excitation. In the initial 
phase of amplified spontaneous emission (pft = ... 1) 
one can see that the radiation spectrum is narrower than 
the inhomogeneous atomic spectrum. This gives the first 
indication of spectral condensation, a well-known phe- 
nomenon in amplifying media and lasers. 



B. Effects of radiative atom-atom interactions on 
intrinsic optical bistability 

One of the most interesting dynamical effects in dense 
media due to the Lorentz-Lorenz nonlinearity is thepos- 
sibility of intrinsic optical bistability predicted in M. If 
a radiatively broadened two-level system is resonantly 
driven by a coherent field of Rabi-frequency fl it shows 
mirrorless, i.e. intrinsic bistability, if the atomic density 
exceeds some critical value. The bistability results from 
an effective feedback introduced by the Lorentz-Lorenz 
correction. 

We here consider a dense ensemble of resonantly driven 
two- level systems as shown in Fig. 9. For simplicity of 
the present discussion we assume that the driving field is 
homogeneous. 




Y* Y 



£2l 



-1 T 1 



FIG. 9. Strongly driven, dense two-level medium. 7 and 7* 
describe radiative and non-radiative decays respectively, and 
Ql denotes Rabi-frequency of local field. 

For most practical realizations this assumption is not 
valid. We are here however interested in principle ques- 
tions and will therefore ignore drive-field depletion. The 
assumption of a homogeneous driving field implies a ho- 
mogeneous behavior of the atomic system and we can 
disregard the spatial dependence in the collective decay 
and light-shift terms. The density matrix equations for 
the system under consideration read in a rotating frame 

Paa = -TaPaa + T 'pbb - i(^Pab - ^Pafc), (83) 
Pab = -^abPab - + Cjp a b){p a a ~ Pbb), (84) 

where we have assumed a real 51. r a = 7 + 7* + T is 
the total population decay rate out of the excited state, 
with 7 and 7* being the free-space radiative and non- 
radiative decay rates, and T the collective decay rate. 
r fl 6 = T + (7 + 7*)/2. There is no collective light-shift 
contribution here due to symmetry reasons. One recog- 
nizes a term proportional to the atomic polarization p a b 
that adds to the Rabi-frequency fl. This term is due to 
the Lorentz-Lorenz correction ( |3l| ) and has the charac- 
ter of a feedback (atomic polarization generates a field 
contribution Cjp a b which couples back to the atom). 
C = NX 3 /4tt 2 is the cooperativity parameter, that es- 
sentially determines the number of atoms in a volume 
A 3 . 

The stationary solution of Eqs.(|8^) and ( |S4| ) for the ex- 
cited state population for T = 7* = is shown in Fig. 10 
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for different cooperativities. One recognizes bistability 
for C> 3. 



P 0, 



C=l. 



Q/y 

FIG. 10. Stationary excited state population as function 
of driving-field Rabi-frequency £1 for different cooperativity 
parameters. Here V — 7* = 0. 

We now analyze the effect of incoherent photons in- 
side the medium. To self-consistently determine the col- 
lective decay rate, we have to calculate correlations of 
the dipole operators a = \b)(a\ and cA in terms of den- 
sity matrix elements. More precisely we need the sec- 
ond order cumulants, i.e. correlations of operators mi- 
nus their mean values (which are nonzero in the present 
case) . Thus we start with the Heisenberg-Langevin equa- 
tions for Scr = a — (a) and Sa aa — a aa — (a aa ), where 
Oaa = \a)(a\: 

6a = —T a b6a — 2iQ,i J 8a aa + noise (85) 
S<J aa = -2T ab Sa aa - i(ft* L Sa - Q L SaA + noise, (86) 

where Ql = f2 + C"/p a b- The relevant correlations 
can be obtained from these equations by Laplace- 
transformation. This yields at resonance (u = uj a b) 



P s (u ab ) 



p 2 N T ab (p aa - p bb ) 



r 2 

L ab 



2in, 



(87) 



2p 2 N 



ab^Paa 



\Pab\ 



2\Vl L \ 2 p aa pbb 



T ab [ri b + 2\n L \ 



Note that we have omitted space and time arguments, 
since we are interested in the stationary properties of the 
system in a homogeneous sample. From ( p7|) we imme- 
diately find the absorption coefficient 



P 2 N UJ ab 



3heo c 



r 2 

L ab 



2 0, 



(Pa 



Pbb) 



(89) 



We now consider a thin plate of thickness d as shown in 
Fig. 9 and assume that the beam diameter of the driving 
field is large compared to d. Carrying out the spatial 
integrations in Eq.(|9|) using fl89|) and ( |8"S| ) we find the 
following relation for the collective decay rate 



r 

7 



(90) 



Pbb — Paa 



\Pab\ 



2\n, 



F 2 

1 ab 



-PaaPbb 



This equation for T is not yet explicit since T enters the 
right hand side of the equation in an essentially nonlin- 
ear way. Using the stationary solutions of the density 
matrix equations (^) and ( |S4|) with T as independent 
variable, one can (with some additional approximations) 
solve Eqs.©. This yields 



F 

7 



Paa 



l-e 



-K 



with 

K = Cr 



1 - 2p aa 

l+ 7 */7(l-2p aa ) 



(91) 



2% + 2c 2 (i - 2 Paa y + {1+ ^CX1T ))2 



(92) 



C is the cooperativity parameter, and r = ird/\. 

In Fig. 11 we show the stationary solutions for the 
excited state population as function of the driving-field 
Rabi-frequency ft for different cooperativity parameters 
and for purely radiative decay, i.e. 7* = 0. The dot- 
ted curves correspond to the solutions without radiative 
atom-atom interactions. As can be seen, bistability per- 
sists, but cannot be resolved for physically reasonable 
values of r (r = 100 in Fig. lib). Radiation trapping 
prevents the energy to escape from the sample and al- 
ready very small external pumping is sufficient to keep 
the atoms in a highly excited state. 




Rm 0.4 




FIG. 11. Stationary excited state population as function of 
driving-field Rabi-frequency fi for radiative decay (7* = 0) 
and with radiative atom- atom interactions, r = nd/\ = 1 
(upper picture) and r = 100 (lower picture). Dotted curves 
show behavior without radiative atom-atom interactions (see 
Fig. 10.) 

The situation is different, if there is also non-radiative 
decay, as shown in Fig. 12 (here 7* =7). In comparison 
to the radiatively-broadened case of Fig. 11, the bista- 
bility curves are only moderately altered even for large 
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samples with r — 1000. The critical cooperativity at 
which bistability starts to occur is somewhat increased. 
The non-radiative decay provides an additional energy 
escape channel, such that the trapped incoherent radia- 
tion is not strong enough to keep the atoms in the excited 
state. 

Thus we can conclude, that radiative atom-atom in- 
teractions do not destroy intrinsic bistability in driven 
two-level systems, if non-radiative decay is present. This 
is different from our previous result 291. The reason for 
this discrepancy is, that our previous approach essentially 
neglected the medium effect on the retarded propagation 
and was therefore inconsistent for larger densities. This 
would correspond to replacingthe retarded propagator 
of the interacting field in Eq.(|45|) and Fig. 3 by the cor- 
responding free-space Greensfunction. 
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FIG. 12. Stationary excited state population as function of 
driving-field Rabi-frequency £1 in the presence of additional 
nonradiative decay (7* = 7) and with radiative atom-atom 
interactions, r — 1000. Dotted curves show corresponding be- 
havior without radiative atom-atom interactions for 7* = 7. 

We note that we did not intend to present a compre- 
hensive discussion of effects that could affect intrinsic op- 
tical bistability. In particular in atomic vapors collisions 
may have a much more pronounced effect. Furthermore 
the depletion of the pump field needs to be taken into 
account and realistic experimental schemes such as selec- 
tive reflection spectroscopy [£7j need to be considered. 

VI. SUMMARY 

We have shown that the interaction of a classical ra- 
diation field with a dense ensemble of atoms can be 
described by modified Maxwell-Bloch equations in the 
Markov- limit. While the equations for the macroscopic 
classical field, the Maxwell equations, remain unchanged, 
the atomic equations of motion have additional nonlin- 
ear and spatially nonlocal terms that result from the ex- 
change of spontaneous photons between the atoms. The 
first nonlinear term, the LL-correction, is only present 
if there is initial coherence or an external coherent field. 
The nonlinear and nonlocal collective decay and level- 
shift terms are present whenever there is excited-state 
population. In the Markov-limit of spectrally broad inco- 
herent radiation inside the medium, the modified Bloch 



equations have the form of single-atom density matrix 
equations. They are of the Lindblad type and thus fulfill 
the formal requirements for conservation of probability 
and positivity. The effect of the collective decay has been 
discussed for the examples of an initially excited system 
of inhomogeneously broadened two-level atoms and in- 
trinsic optical bistability. In the first case we find accel- 
erated decay (superluminescence) in the initial phase and 
radiation trapping in the final phase. In the long-time 
limit, where the excited state population is small, the 
equations of motion coincide with the Holstein equations 
of radiation trapping Jl3| . The collective decay modi- 
fies intrinsic optical bistability in a strongly driven two- 
level system. As opposed to our previous prediction Q , 
bistability persists however, if also non-radiative decay is 
present. Collective decay and pump processes as well as 
light shifts are relevant for the population dynamics and 
are particularly important for ground state coherences in 
multi-level systems. A detailed discussion of coherence 
effects in dense multilevel systems as well as the study of 
cooperative decay processes will be the subject of future 
work. 



ACKNOWLEDGMENT 

The authors would like to thank G. Agarwal, 
C. M. Bowdcn, J. P. Dowling, M. Lukin, A. Manka, 
A. Schenzle and M. O. Scully for stimulating discussions. 
The support of the Office of Naval Research and the 
Texas Advanced Research and Technology program are 
gratefully acknowledged. M.F. would also like to thank 
the Alexander-von-Humboldt foundation and S.F.Y. the 
Studienstiftung des Deutschen Volkes and the German 
Academic Exchange Service (through their program HSP 
II) for financial support. 

Appendix A: Free-space decay rate 

In Eq. (|2l]) we have given the spontaneous decay rate in 
the atomic medium. We now show that this expression 
leads to the well-known Wigner-Weisskopf result for the 
radiative decay of a two-level atom if the interacting field 
is replaced by the free field. For a single two-level transi- 
tion with dipole moment along we have according to 

(1) 

7=^ j M jT([E^(f ,t + T),E^(f ,t)})e^. (93) 

The free-field commutator is given by |22| 

[S + M (f ,t + r),£ - M (fo,t)] - (94) 
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Substituting (||) into @ yields 

P 2 ^ 1 f 2 / fc 2 \ 

Finally carrying out the angle-integration leads to 

SttV 



with 



7 



p 2 w 3 



37r?ieoc 3 3?ieoA 3 ' 



(95) 



(96) 



n#(i, 2) = n++(i, 2) - n+-(i, 2) = 
n-+(i,2)-n;-(i,2) = 



(104) 



which is the free-space spontaneous emission rate from 
the Wigner- Weisskopf theory [E5J . 



Appendix B: Dyson equation for real-time 
Greensfunctions 

In this Appendix we derive the integral equations (|4^) 
and ( |46"| ) for the real-time GF from the Dyson equation 
(||) for the contour GF. Noting that L>~+ = in RWA, 
we immediately find from (pt2|) 



Thus we obtain the Dyson-equation (46) for the retarded 
Propagator inside the medium of Sec. IV. 

We now turn to D~ + . Substituting II — = rr+-IT ot 
in (pOh we find 



D 



(i 





n+ 


+ D++- 


f-£> + 


+ IT 1 


- D- 




n 


-D-+- 




n 


+ D+ 




rr 


-D-+- 






+ D+ 



D-+ = 

where we used a short notation DTs 



(97) 
(98) 



II— D~ 



Jfd2dAD^(l,3)U--(3A)D-+(i,2) and the Integra- 
tion goes over physical times from — oo to oo and over the 
volume of the sample. Note the sign changes in Eqs.(|9^) 
and ( |9~7j ) resulting from the fact that the contour inte- 
gration on the lower branch goes in the reverse direction. 
Making use of d3S|-[4l|) we obtain 



D-+ = -D 1 ct n rct D-+ - £>{, ct ir+L> adv . (105) 
Iteration of this equations yields 

D~+ = - i>;; : - DinVl);:' (106) 

ir+L> adv 

which can be rewritten in the compact form 

D-+(l,2) = -J J d3d4^(l, 3) n^(3, 4)^(4,2), 

(107) 



^rctjrrct^rotjrrot^rot 



where 



7-jadv 



IT 



nadv 



D 



-+ 



IF 



-ir+Jz^ 

H-+D++- 



u' (i,2) = n-+(i,2) 



(108) 



E((<4(*i)M*2))) 6 ft - ? i) 5 ft - f i) 



D-+ = -DS 



n 



-+ n+ 



- ir 



D~ 



, (99) 
(100) 



Thus we arrive at Eq.(|45|) of Sec. IV. 



Applying the relation 

n adv (i, 2) = n++(i, 2) - ir+(i, 2) = 
n +- (i,2) - rr~(i,2) 



(101) 



and subtracting Eqs.(|100j) from ( p9|) yields the Dyson 
equation for £> adv = D++ - £>-+: 

B#(l, 2) = D%%(1,2)- (102) 

J J d3dADf^(l,3)nfj(3,A)Dlf(4,2) 

where we have restored full notation. Since Dj^(l,2) = 
D adv (2,l) one immediately obtains the corresponding 
Dyson equation for the retarded propagator 

D%(1,2) =DI%(1,2)- (103) 
J J MdAD&^l, 3)^^(3,4)^(4, 2) 



Appendix C: Solution of the Dyson equation for D TCt 

With the approximations made in Sec.IV-3 we derived 
the matrix solution ([55]) for the Fourier-transform of the 
retarded GF in the medium. 



D {q,u,t) 



1+ D ( 9 »-P ret (o;,t) 



- 1 —ret 

•D (q,u). 
(109) 



To evaluate this expression we first approximate the 
Fourier-transform of the free-space retarded propagator. 
According to p3] 



(110) 



if i 



9(r) 



6 92 c 2 92 



5(r — cr) 



13 



with r 



\fi — r*2| and t = t± — £2- Thus 

ifi fu 2 . x a X/3 d 2 \ e~ lujr / c 



Aire 



r 2 dr 2 



(111) 



For large u, such that A <C r, only the spatial derivative 
of the exponential contributes and we find 



iTiuj 2 ( 



-iujr I ' c 



(112) 



We now approximate ( |112| ) by ignoring the polariza- 
tion, i.e. by performing an orientation average. 



Xa.Xf3 



1 



Oa/3- 



(113) 



This approximation is exact when the medium is ran- 
domly polarized. (113) leads to 



ur/c 



6neoc 2 



and thus 



D {q,u) = - 



3e c 2 q 2 - 4 + 2it u - 



(114) 



(115) 



where we have introduced a small positive constant e to 
move the pole at q = lu/c into the lower half of the com- 
plex plane. 



Substituting ( |115| ) into ( |109| ) we find 



D rct (x,iu) 



,. , d?q D (Me-*** 
Zir ) J 



(116) 



hio 2 



6ir 2 eoc 2 r 



dq 



qe 



1 



2ih 
3eo 



P ret (w) 



2ie L - 



We evaluate this integral by contour integration along 
the real g-axis and back in the lower half plane (note that 
r > 0). In free space (P = 0) and for e = we would 
have two poles on the real axis at q = ±lu/c. In the 
free-space case we had to introduce the constant e > 
to move the pole at q — lu/c into the lower half-plane 
(and at the same time to move the pole at q = —lu/c 
into the upper one). This is necessary to have retarded 
propagation. For e < one would obtain the advanced 
propagator. If an absorbing medium is present, Re[P ret ] 
is negative and hence adds to e. Therefore again the 
only contributing pole is that at q w ui/c. However, if 
the medium is amplifying there is a problem. In this case 
Re[P ret ] is positive and counteracts e. In such a case the 
pole at q = —lu/c could move into the lower half-plane 
and we would obtain an advanced instead of a retarded 



propagator. The origin of this problem is, that the re- 
tarded propagator in an amplifying medium is strictly 
speaking not Fourier-transformable, since it is an expo- 
nentially growing function of distance. In such a case one 
has to take into account the finite spatial dimensions of 
the amplifying medium and introduce a cut-off function 
which leads to a finite-space Fourier-transform. Although 
this lacks mathematical rigor, we now assume that the 
effect of the cut-off function is modeled by a sufficiently 
large value of e, such that the pole at q as lu/c remains in 
the lower half-plane. With this we find 



D Tet (x,Lu) = 



where 
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